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The critical behavior of the two-dimensional O(N) model close to criticality is shown to be en- 
coded in the fractal structure of the high-temperature graphs of the model. Based on Monte Carlo 
simulations and with the help of percolation theory, De Gennes' results for polymer rings, corre- 
sponding to the limit TV — > 0, are generalized to random loops for arbitrary —2 < N < 2. The loops 
are studied also close to their tricritical point, known as the O point in the context of polymers, 
where they collapse. The corresponding fractal dimensions are argued to be in one-to-one corre- 
spondence with those at the critical point, leading to an analytic prediction for the magnetic scaling 
dimension at the O(N) tricritical point. 



The high-temperature (HT) representation of the crit- 
ical 0(A) spin model naturally defines a loop gas, cor- 
responding to a diagrammatic expansion of the partition 
function in terms of closed graphs along the bonds on 
the underlying lattice 0. In the limit N — > 0, the loops 
reduce to closed self-avoiding random walks first consid- 
ered by de Gennes as a model for polymer rings in good 
solvents at sufficiently high temperatures, so that the van 
der Waals attraction between monomers is irrelevant . 
In his seminal paper, de Gennes related the fractal struc- 
ture of self-avoiding random walks to the critical expo- 
nents of the 0(N — > 0) model. Invoking concepts from 
ercolation theory and recent Monte Carlo (MC) data 
U of the HT representation of the two-dimensional 
(2D) Ising model (N = 1), we extend in this Letter de 
Gennes' result to arbitrary —2 < N < 2 in 2D. We con- 
sider the loops also close to their tricritical point where 
they collapse, known as the © point in the context of 
polymers. We argue that the 2D fractal dimensions at 
the tricritical point are in one-to-one correspondence with 
those at the critical point, allowing us to also predict the 
magnetic scaling dimension at the 0(A) tricritical point. 
We support our theoretical prediction by comparing it 
with recent high-precision MC data |5j . 

A particularly simple representation of the O(A) uni- 
versality class is specified by the partition function 
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where the product is over all nearest neighbor pairs, and 
the spins have N components S x = . . . , S£ ) and 

are of fixed length S x = N. The trace Tr stands for 
the sum or integral over all possible spin configurations. 
The weighting factor is obtained by truncating the more 
standard Boltzmann weight exp(AS x • S x ')- This choice 
mimics the weighting factor of the Ising model, where 
exposes*/) oc (1 + KSxStf) with K = tanh/3. When 
formulated on a honeycomb lattice, which has coordina- 
tion number z = 3, the HT graphs of the truncated model 
are automatically nonintersecting and self-avoiding. The 
partition function can then be written simply as a sum 



over all possible closed graphs p, Z = J2{g\ K h N l , with 
b and / the number of occupied bonds and separate loops 
forming the graph G. The parameter K in the spin for- 
mulation Q appears as bond fugacity in the loop model. 
By mapping it onto a solid-on-solid model, the critical 
exponents as well as the critical point were determined 
exactly 0. 

In the high-temperature phase, the HT graphs have a 
finite line tension 9q and are exponentially suppressed. 
A typical graph configuration in this phase shows only a 
few small loops scattered around the lattice. Upon ap- 
proaching the critical point from above, the lattice starts 
to fill up with more and also larger graphs. At the crit- 
ical point, the line tension vanishes, causing the expo- 
nential suppression to disappear. Graphs of all sizes now 
appear in the system as they can grow without energy 
cost, i.e., the HT graphs proliferate. A graph spanning 
the lattice can be found irrespective of the lattice size- 
much like the appearance of a spanning cluster at the 
percolation threshold in percolation phenomena Q . The 
average number density £b of graphs containing b bonds 
takes asymptotically a form similar to that of clusters in 
percolation theory, 
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with gq and tq two exponents whose values define the 
universality class. The line tension vanishes upon ap- 
proaching the critical point at a pace determined by the 
exponent ctq- When present, this Boltzmann factor ex- 
ponentially suppresses large graphs. The algebraic factor 
in the graph distribution is an entropy factor, giving a 
measure of the number of ways a graph of size b can be 
embedded in the lattice. The configurational entropy is 
characterized by the exponent tq. As in percolation the- 
ory |3, it is related to the fractal dimension Dq of the 
HT graphs via 



r G = d/ Dq + 1, 



(3) 



with d — 2 the dimension of the lattice. 

When summed over all sizes, the graph distribution 
yields the scaling part of the logarithm of the partition 



2 



Model 


N K- 


c 


a (3 j r] v 


A3 TQ 0"G 


Gaussian 
SAW 
Ising 
XY 


-2 \ 

§ 

1 ! 

2 1 


-2 


i 

2 
1 


1 1 | 

1 5 43 5 3 

2 64 32 24 4 

i - - 1 

u 8 4 4 - 1 - 

— OO OO OO j OO 


5 13 8 
4 5 5 
4 5 1 
3 2 

11 27 8 
8 11 11 
3 2 q 
2 3 u 



TABLE I: Critical exponents of the two-dimensional criti- 
cal O(N) spin model, with N = —2,0, 1,2, respectively, to- 
gether with the parameter k_ and the central charge c = 
1 — 6(1 — k_) 2 /k_. In addition, the fractal dimension Dq of 
the HT graphs as well as the two exponents tg, o"g character- 
izing their distribution are given. For N = 1, the values of the 
latter two were recently established numerically in Ref. Q. 

function, 

InZ~^4- (4) 

b 

Each graph therefore contributes equally to the scaling 
part of the free energy, irrespective of its size. 

Table [I] summarizes the critical exponents and fractal 
dimensions of the four most common O(-ZV) models. The 
negative value N — — 2 corresponds to the noninteract- 
ing model, for which the critical exponents take their 
Gaussian values 0. In the polymer limit N — > 0, first 
studied by de Gennes 0, the fractal dimension of the 
HT graphs is simply the inverse of the correlation length 
exponent v. In general, however, it follows from Tabled 
that D G ^\jv. To generalize de Gennes' result, we note 
that in percolation theory a similar relation between the 
fractal dimension D of clusters and the correlation length 
exponent involves the Fisher exponent a, viz. 

D = 1/crv. (5) 

A closer look at de Gennes' derivation reveals that o"g = 1 
in that case, implying that the result for polymers in good 
solvents or self-avoiding random walks is consistent with 
Eq. (0) . In a recent MC study of the HT graphs of the 
2D Ising model , we numerically found the value ctg = 
0.732(6). This estimate is within one standard deviation 
from the value ctq = 8/11 = 0.7273... expected from 
Eq. (JSJ , with v = 1 and the fractal dimension Dq = 11/8 
appropriate for the Ising model. 

Parameterizing the 2D 0(N) models as [Tol [TH iV = 
— 2 cos(7r/S„), with | < re_ < 1, we obtain from Eq. (JSJ) 
ctq = 8(1 — k_)/(2 + R-), where use is made of the known 
results l/u = 4(1 - /t_) and [13 

D G = l + S_/2. (6) 

The entropy exponent © follows similarly as tq — 
(6 + k_)/(2 + yielding tq = 5/2 for a self-avoiding 
random walk and tq — 27/11 = 2.4546 . . . for the Ising 
model. Through the exact enumeration and analysis of 



the number of self-avoiding loops on a square lattice up 
to length 110, the expected value tq = 5/2 has been 
established numerically to very high precision I n 
our MC study 3] , we numerically obtained the estimate 
tq = 2.458(5) for the Ising model in good agreement with 
the theoretical prediction. 

The O(N) spin-spin correlation function G(x, x') = 
G(x — x') is represented diagrammatically by a modified 
partition function, obtained by requiring that the two 
sites x and x' are connected by an open HT graph T . 
On a honeycomb lattice, the scaling part of the correla- 
tion function is given by the connected graphs 

G(x, x') ~£^ = E *6(x, ^)K\ (7) 
{r} b 

where Zf,(x, x') is the number of (open) nonintersecting 
and self-avoiding graphs along b bonds connecting x and 
x'. It is related to the graph distribution J5J through 
lb = (l/^)Ex z ' | ( x ' x )' w ^ u ^ the lattice volume. 
Since Zb(x, x) refers to closed graphs starting and ending 
at x, the factor 1/b is included to prevent overcounting as 
a given loop can be traced out starting at any lattice point 
along that loop. Strictly speaking, G(x, x') < J2{r\ K b 
as the cancellation of the disconnected graphs in the nu- 
merator and Z in the denominator, required for an equal- 
ity, is not complete: For a given open graph T, certain 
loops present in Z are forbidden in the modified parti- 
tion function as they would intersect T, or occupy bonds 
belonging to it, which is not allowed. In other words, the 
presence of an open graph influences the loop gas and vice 
versa. Since each loop carries a factor N, the loop gas 
is absent in the limit N — > 0. The inequality in Eq. (JJJ 
then becomes an equality and the open graphs become 
ordinary self-avoiding random walks with Dq = 4/3. For 
N > 0, the loops obstruct the formation of graphs con- 
necting the two endpoints, so that the fractal dimension 
of these self-avoiding graphs on the honeycomb lattice is 
larger. 

The magnetic susceptibility x follows as \ = 
£ x ,G(x,x') ~ E b ZbK\ with z b = £ x ,^(x,x') the 
number of open graphs of size b starting at x and ending 
at an arbitrary lattice point. For the susceptibility to di- 
verge with the correct exponent 7, it must behave close 
to the critical point as 

b 

for large b, where like in the closed graph distribution J2J, 
the Boltzmann factor suppresses large graphs as long as 
the line tension #g is finite. Indeed, replacing the sum- 
mation over the HT graph size b with an integration, we 
find x ~ \K— . The asymptotic form (JSJ generalizes 
the de Gennes result for N — * with ctq = 1 to arbitrary 
—2 < N < 2 with (Tq / 1. Note that on account of 
Fisher's scaling relation 7 = (2 — rf)v, the combination 
ctg7 m Eq. (JSJ satisfies ctg7 = (2 — T])/Dq. 
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The ratio of z&(x, x') and z& defines the probability 
Ph(x, x') of finding a graph connecting x and x' along b 
bonds yjj. On general grounds, it scales at criticality as 
(d=2): 



/£>G 



F b (x, x') = z 6 (x, x')M ~ b- d ' DG P (|x - x'l/fe 1 

'(9) 

with P a scaling function. Since at the critical point, 
ZhK\ ~ ^gt-i according to Eq. (JSJ), we obtain for the 
correlation function 

G(x,x') ~ 5> b P b (x,x')A c b ~ l/|x-x'|^ 2+ ", (10) 
b 

where use is made of Eq. iJSJ and Fisher's scaling relation. 
Equation (|10|l is the standard definition of the critical 
exponent rj, whose exact value is given by 



r) = 2- 3/4k_ - k_, 



(11) 



and thus provides a consistency check. Also, with v = 
(t g - l)/dcr G , as is implied by Eqs. © and ©, Eq. (gj 
yields the scaling relation dv = 2 — a, where a deter- 
mines the scaling behavior of the free energy close to the 
critical point, InZ ~ |if — A c | 2 ~ a . Apart from v and a, 
the exponents have no simple dependence on the graph 
distribution exponents (Tq, t g- This is because the oper- 
ator whose scaling dimension tjq is given by the fractal 
dimension Dq of the HT graphs is not a simple one, con- 
sisting of two spin components at the same site S^S^ 
which measures the tendency of spins to align • 

In our MC study Q , rather than determining the scal- 
ing dimension y G directly, we measured the so-called per- 
colation strength P^ , giving the fraction of bonds in the 
largest graph. This observable obeys the finite-size scal- 
ing relation ~ L~^ G ^ V , where the exponent (3q is re- 
lated to yc via y G — d — @g/v. We found (3 G — 0.626(7) 
for the Ising model in perfect agreement with the value 
(3q = 5/8 = 0.625, leading to y G — 11/8, which coincides 
with the fractal dimension Dq of the HT graphs. 

We next extend our results for the HT graphs at the 
critical point to the tricritical point. It is generally ac- 
cepted that including vacancies in the O(N) model gives 
in addition to critical behavior rise to also tricritical be- 
havior. By gradually increasing the activity of the vacan- 
cies, the continuous 0(A) phase transition is eventually 
driven first order. The endpoint where this happens is 
a tricritical point. In the context of polymers (A — » 0), 
the latter obtains by lowering the temperature to the 
point where the increasingly important van der Waals 
attraction between monomers causes the polymer chain 
to collapse. Coniglio et al. 0] argued that a polymer 
ring at the O tricritical point is equivalent to the hull 
of a percolation cluster. With the known dimension for 
the percolation hull 0| , the analogy then gives D G = | 
as fractal dimension for a polymer chain at the point 
(the superscript "t" refers to the tricritical point) . It also 



implies that polymer chains at the critical and tricritical 
point share the same central charge c because both the 
0(N — > 0) model and percolation have c = 

To generalize the analogy found by Coniglio et al. ^5| , 
we consider the Q-state Potts model, which in the limit 
Q — » 1 describes ordinary, uncorrelated percolation 
(lif. In the Fortuin-Kasteleyn (FK) formulation [l9j . 
the model with Q > 1 is mapped onto a correlated per- 
colation problem. Clusters are formed by lumping to- 
gether with a certain temperature-dependent probability 
nearest neighbor spins in the same spin state. These so- 
called FK clusters percolate at the critical temperature 
and their fractal structure encodes the entire critical be- 
havior. In 2D (and in 2D only), also the geometrical clus- 
ters, formed by unconditionally lumping together nearest 
neighbor spins in the same spin state, percolate at the 
critical temperature. Their fractal structure encodes the 
tricritical Qt-state Potts behavior, which emerges when 
enlarging the pure model to include vacancies. The tri- 
critical point shares the same central charge c with the 
critical point, but apart from c = 1, Q(c) ^ Qt(c). Both 
fractal structures and thus both critical behaviors are in- 
timately related, being connected by replacing k_ with 
k+ = 1/R- > 1 in the a ppr opriate expressions (for de- 
tails, see Ref. Q or Ref. |2(j. where similar conclusions 
were reached independently). This map conserves the 
central charge. With increasing Q or c, the critical and 
tricritical points approach each other until merging at 
c = 1, where Q = Qt — 4 and the scaling behaviors of 
FK and geometrical clusters coincide. Beyond Q = 4, 
the transition becomes discontinuous. 

The hulls of the geometrical clusters of the Q(c)-state 
Potts model at the same time represent the HT graphs of 
the critical O(N) modeh with < N(c) < 2, sharing the 
same central charge |3, |2lj . Since the geometrical clus- 
ters encode the tricritical Potts behavior, while the FK 
clusters encode the critical Potts behavior (characterized 
by the same central charge), it is natural to expect the 
FK hulls to represent the HT loop gas not at the criti- 
cal, but at the tricritical point. For the special case of 
polymers (N — > 0), this reproduces the result by Coniglio 
et al. 15] . Note that in both Potts and O(A) models, 
including vacancies leads to tricritical behavior. How- 
ever, the critical and tricritical points get interchanged 
when passing from one model to the other. The fractal 
dimension D G of the tricritical loops follows by applying 
the central-charge conserving map k_ — > R + — 1/R- to 
Eq. ©. We submit that the exponent rjt follows in the 
same way from Eq. Given rj t , the scaling relations 

then yield values for the ratios Pt/^t and 7t/^t- 

To verify these predictions we rewrite our results, given 
as a function of k_, as entries in the Kac table 



[(to + l)p — mq] 



1 



(12) 



p - q 4to(to+1) 
Here, the parameter m is related to k_ and the central 
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FIG. 1: Our analytic prediction for the tricritical exponent 
r/t (solid line) as a function of the central charge c compared 
with high-precision MC data (marks with error bars) 0. Our 
prediction extends only to c = 1, indicated by the vertical 
line, while the numerical results continue to the tricritical 
0(2) model with c = 1.149(1) @. 

charge via k_ = m/(l + m) and c = 1 — 6/m(m + 1), 
respectively, while the central-charge conserving map be- 
comes to — > — m — 1. Specifically, for < c < 1 

D^ = l + l/2/«_=2-2/i m , m , (13) 

leading to the correct result -Dq — 2 — 2/12,2 = | for 
a polymer at the point (m = 2), and 

?7 t = 2 - 3k_/4 - 1/R_ = 4/i m/2!m/2 . (14) 

The right hand of this equation is in agreement with the 
results T] t = 4/11,1 = for the 9 point (to = 2) 
and rjt = 4/i 2 ,2 = 55 f° r the tricritical Ising model 
(m = 4). After circulating a draft of this paper, we have 
been informed about a recent high-precision MC study of 
the tricritical O(N) model in Ref. [5j previously unavail- 
able to us, in which the authors extend earlier numerical 
work on the tricritical 0(|) model 23] . In Fig. ^ we 
compare our theoretical prediction (|14fl for 7y t with the 
magnetic scaling dimension x\ = 2 — y\ = i] t /2 obtained 
numerically in that study. For < c < jq, the MC data 
are within one standard deviation of our prediction. Be- 
yond the tricritical Ising model (c = jq) the numerical 
data start deviating from our analytic result. A detailed 
future investigation is required to clarify this discrepancy. 

The tricritical HT graphs, representing simultaneously 
the hulls of FK clusters, have a distribution again of the 
form (J2J, characterized by two exponents CTq,Tq. Given 
our result 1|13[) for the fractal dimension of these graphs, 
Tq is determined exactly through Eq. © with tq and 
Dq replaced by their tricritical counterparts Tq and D l G . 

In conclusion, we have shown that for — 2 < N < 2 the 
fractal structure of 2D HT graphs of the 0{N) spin model 
encodes the 0{N) critical behavior. We thereby extended 
de Gennes' result for self-avoiding loops in the limit 
N — > to random loops for arbitrary —2 < N < 2. We 



studied the loops also close to the point where they col- 
lapse, corresponding to the HT representation of the tri- 
critical O(N) model. The fractal structure of the tricrit- 
ical loops was argued to be in one-to-one correspondence 
with that of the critical loops, allowing us to also predict 
the magnetic scaling dimension x\ = %/2 = 2/i m / 2 , m /2 
at the O(A^) tricritical point, in very good agreement 
with recent MC data in the range < c < 0.7. 
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